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The quarks in quenched QCD have an anomalous self-interaction in the flavor singlet Goldstone 
boson channel. This coupling is extracted from a graph with disconnected quark lines, and is used to 
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action, the Witten-Veneziano relation is an exact relation between the topological susceptibility (as 
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action we compute the hairpin amplitude and determine the fermion zero-mode susceptibility, the 
inferred eta-prime mass and other parameters characterizing the low energy chiral properties of 
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I. INTRODUCTION 



The physics of the flavor singlet pseudoscalar eta-prime meson has long been a source of puzzles. 

• The eta-prime is not a Goldstone boson because of the presence of the axial anomaly. 

• The eta-prime propagator in full QCD involves a set of connected and disconnected quark-line diagrams. Sum- 
ming a plausible subset of these diagrams (see Fig. [l]) produces a shift in the eta-prime mass away from the 
masses of the true Goldstone bosons. But what gauge field dynamics "fills in the white space" between the 
quark loops? 

• In the limit of a large number of colors N c , the eta-prime mass is expected to scale as 1/N C as opposed to the 
masses of non-Goldstone mesons, which remain order one. A mass formula relating the eta-prime mass to the 
topological susceptibility has been derived by Witten and Veneziano JlJ 

m\, + ml - 2m 2 K = $ = ^ (1) 



where \t is the topological susceptibility of the pure gauge theory, Nf is the number of flavors, and the pion 
decay constant, defined through (O|-07o75V'l 7r ) = mirfwi has the experimental value of 132 MeV. The right-hand 
side is formally 0(1/N C ) because / w ~ \fW c - The suppression of fermion loops at large N c is similar to the 
absence of fermion loops in the quenched approximation, and lattice calculations of the topological susceptibility 
in quenched QCD (using pure gauge observables) give numbers with good numerical agreement with Eq. (Q), 
when evaluated with the physical masses of the particles. It seems strange that this large- N c formula should 
be quantitatively correct, since the eta-prime mass is not particularly small compared to the masses of other 
non-Goldstone boson mesons. (It is worth remarking that the flavor-nonsinglet pseudoscalar sector is also the 
home of other large-mass quantities: one example is the scale of chiral symmetry breaking: m^/(m u + rrid) — 3 
GeV. The physical eta-prime mass is less than one third of that number.) 

The relation of these statements to each other is an ongoing source of controversy. Lattice simulations can in principle 
contribute to the discussion and indeed, there have been many studies of the eta-prime on the lattice, dating back 
to the earliest days of simulation (a partial list includes [^]-||). Most of these studies have been in the quenched 
approximation, where the mass of the eta-prime is inferred from the size of the hairpin graph compared to the 
connected graph, and compared to the topological susceptibility. When the contribution of individual eigenmodes of 
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the Dirac operator to the hairpin graph is computed, it is often seen that the low-lying eigenmodes make a substantial 
contribution to the hairpin graph. 

This paper is another calculation of the hairpin diagram in quenched QCD. Why revisit this question yet again? 
The reason is that the calculations presented here are done with a fermion action which respects chiral symmetry on 
the lattice via the Ginsparg- Wilson relation |)|: an overlap JlQ] action. All previous studies of the eta-prime were 
done with lattice actions which have chiral symmetry breaking artifacts: for Wilson-type fermions, for example, real 
eigenmodes of the Dirac operator are not zero modes nor are they eigenstates of 75. Since the properties of the 
eta-prime are known to be intertwined with the anomaly, one would expect that calculations with exact lattice chiral 
symmetry might be more revealing for the study of the eta-prime. 

Indeed that is the case: when the fermions are described by an overlap action, the Witten-Veneziano relation is an 
exact relation between the expectation value of the hairpin diagram-which is given by the susceptibility of fermion 
zero modes-and the anomalous coupling between two flavor singlet Goldstone bosons. Interpreting that coupling as 
the eta- prime mass and equating the zero mode susceptibility to the topological susceptibility yields Eq. ([!]). 

To see this, first recall some facts about the overlap action: 

The eigenmodes of any massless overlap operator D(0) are located on a circle in the complex plane of radius xq 
with a center at the point (xo,0). The corresponding eigenfunctions are either chiral (for the eigenmodes with real 
eigenvalues located at A = or A = 2xq) or nonchiral and paired; the two eigenvalues of the nonchiral modes are 
complex conjugates. The massive overlap Dirac operator is conventionally defined to be 

D(m q ) = (l-^-)D(0)+m q (2) 

and it is also conventional to define the propagator so that the chiral modes at A = 2xq are projected out, 

D-\m q ) = - l ——{D-\m q ) - -U (3) 

l-m q /(2x ) 2x 

For a summary of useful formulas, see Ref. (for the special case xq — 1/2). 

Now consider the hairpin diagram involving a single flavor of quarks, where the source and sink (black dots in Fig. 
|l|) are the local pseudoscalar density, "tp^ijj. The hairpin diagram then is just 

H(x,y) = (Tv^Dix^y'TT^Diy^y)- 1 ). (4) 

Because only zero modes are chiral, the volume integral of the hairpin graph is proportional to the zero mode 
susceptibility 



where Q is just the number of zero modes, the difference of positive and negative chirality zero modes, Q = n + — n_. 
Regardless of any dynamical model used to describe the hairpin graph, we expect to see a large contribution from 
zero modes to H(x, y) itself, since only they contribute to the susceptibility. 

In quenched QCD, as described by quenched chiral perturbation theory frL2|| , there is an anomalous coupling of two 
Goldstone bosons in the flavor singlet channel, parameterized by a coupling with the dimensions of a squared mass. 
The hairpin graph is analyzed as if each of its quark loops is a propagator for an ordinary pseudoscalar Goldstone 
meson. That is, the momentum space amplitude for the connected graph is 

C(q) = fp^-^fp (6) 
while the hairpin amplitude involving a single flavor is 

q A + m;. Nf <T + 7774 

The quantity /Xq is the squared mass of the "quenched approximation eta-prime" in the chiral limit. (The factor l/Nf 
converts the single-flavor graph into the expectation of the eta-prime mass in iV/-flavor QCD, since each closed loop 
has a multiplicity of Nf, and the wave function (vertex) is scaled by a factor of 1/ ^jNf.) In full QCD the correlator 
which gives the mass of the isosinglet meson is the difference C(t) — NfHf u u(t), and H(t) is supposed to represent 
the first term in a geometric series, which sums up to 
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C(q) - N f H full (q) = C(q) - N f H{q) + . . . = f f 



14^ 



shifting the squared mass of the meson from to m 2 + /ip. In these expressions, fp = (0\ipj5ip\ir) 
from the PCAC relation. Computing the susceptibility directly from Eq. (]?]) gives 

mi N 
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H(x,y) 
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AN f m\ ' 



(8) 

m 2 n U/(2m q ) 



(9) 



and 



we obtain the Witten-Veneziano relation fi^ — 4Nfx/f%, where x is the zero mode 



Equating Eqs. 
susceptibility. 

This simple derivation of the Witten-Veneziano relation from the overlap action uses only the quenched approxi- 
mation, without any reference to the large N c limit. The crucial ingredient is the fact that the hairpin graph, in the 
quenched approximation, takes the form of Eq. (^). Conventional derivations of the Witten-Veneziano relation in 
unqucnchcd QCD proceed rather differently: One begins by considering a correlator of the local topological charge 
density, 



U(k) = / d 4 xexp(ikx)(( 
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16ir 2 N c 



FF(x))( 
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16tt 2 V ( 



■FF(0))) 



which is assumed to be dominated by an eta-prime resonance plus other massive states 

U(k) 



r 2 

k 2 + 77; 



(10) 



(11) 



with C„ 



,,■ - {0\(g 2 /(lQ-K 2 N c )FF\-q'). Using the anomaly equation 2Nf(g 2 /{IQ-k 2 N c )FF = 9 p J 5 M allows one to replace 
the gluonic matrix element by the matrix element of the divergence of the quark current. Computing the susceptibility 
by taking k to zero and setting the eta-prime and pion decay constants equal, so (0|9 M J5W) = \f^l vr ^ ] l fir > gives the 
Witten-Veneziano relation. The difference is that here the eta-prime really is a propagating particle, whereas in the 
quenched approximation the mass is just the value of the two-boson coupling. Note also that in the derivation we 
have just given, the topological susceptibility is replaced by the squared fluctuation in the number of fermionic zero 
modes. The k — ► limit of U(k) is a problematic one due to the contribution of contact terms |l3| ]. 

A third kind of derivation of the Witten-Veneziano relation has recently been given by Giusti, Rossi, Testa and 
Veneziano Jl4[ . Building on derivations of the Ward identity for the flavor singlet axial vector current by Hasenfratz, 
Laliena and Niedermayer and Liischer Jl6|], they show that the Witten-Veneziano relation for overlap fermions 
involves the zero mode susceptibility, not the topological susceptibility. However, they assume that the axial current 
correlator is saturated by a propagating eta-prime resonance, keeping Eq. (|l l[ ) . That does not occur in quenched 
approximation, where the flavor singlet channel propagator is given by Eq. (|7]). 

In full QCD the hairpin correlator is a single particle propagator, but the correlator which gives the eta-prime mass 
is given by C(q) — NfHf u u(q). Hf u u(q) still couples to zero modes, but in order for the difference to couple only to 
the eta prime, the hairpin must have a piece which cancels the connected correlator C(q), the susceptibility of which 
is order 1/m 2 . Then Hf u u(0) — l/m 2 + . . . = x/ m q an d one finds [flT||l7| l the expected result [Q that x — m q- The 
eta-prime mass is only connected to the zero-mode susceptibility through terms which are higher order in the quark 
mass. 
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FIG. 1. A plausible set of quark line graphs, which sum into a geometric series to shift the eta-prime mass away from the 
mass of the flavor nonsinglet pseudoscalar mesons. In the quenched approximation, only the first two terms in the series survive 
as the "direct" and "hairpin" graphs. 

We will compute the zero mode susceptibility both directly and via a calculation of the hairpin graph, and compare 
the results. In practice, in common with all lattice calculations of matrix elements, we will extract ^ from correlators 
using extended sources and sinks. We infer that the mass relation persists, regardless of the choice of source and 
sink, by assuming that the hairpin graph is a two-point vertex of pseudoscalar mesons, as required by quenched chiral 
perturbation theory. 

We will see that Eq. ([5]) produces a measurement of /j, 2 , and % with a somewhat larger statistical uncertainty than 
one would get from a direct measurement (counting zero modes) of (Q 2 )- This happens despite the fact that in the 
fit to the hairpin, one is using information from many points of the lattice, while counting zero modes gives just one 
number per lattice. The signal from the many points just shows very strong correlations from point to point, and the 
actual hairpin graph is computed at nonzero quark mass and must be extrapolated to the chiral limit. 



II. THE LATTICE CALCULATION 



The overlap action used in these studies [ |l9[ ] is built from an action with nearest and next-nearest neighbor couplings, 
and APE-blocked links ^0|. Eigenmodes of the massless overlap Dirac operator D(0) are constructed from eigenmodes 
of the Hermitian Dirac operator H(0) — 75_D(0), using an adaptation of a Conjugate Gradient algorithm of Bunk 
et. al. and Kalkreuter and Simma M]. These eigenmodes are used to precondition the calculation of the quark 
propagator and to construct quark propagators truncated to some number of low lying eigenmodes. 

The data set is generated in the quenched approximation using the Wilson gauge action at a coupling j3 = 5.9. 
The nominal lattice spacing is a = 0.13 fm from the measured rho mass (a value we prefer; see below) or inferred to 
be 0.11 fm from the Sommer parameter using the interpolation formula of Ref. p2| ]. We worked with lattices with 
12 3 x 24 sites. Quark masses in lattice units are am q = 0.02, 0.04, and 0.06, corresponding to pseudoscalar-to-vector 
meson mass ratios of mps/mv — 0.5, 0.61 and 0.67. The fermions have periodic boundary conditions in the spatial 
directions and anti-periodic temporal boundary conditions. We calculated the twenty smallest eigenvalue modes of 
H 2 (Q) in the chiral sector of the minimum eigenvalue, and reconstructed the degenerate opposite chirality eigenstate of 
H 2 (Q) for each nonzero eigenvalue mode. These modes are then recoupled into eigenmodes of D(0). Their eigenvalues 
have imaginary parts ranging up to 0.3/a-0.35/a, or about 500 MeV ]23| ]. 

All correlators we measured will include a smearing function which averages the source of the propagator over a 
localized spatial volume, in order that the meson source resembles an actual hadron wave function. The hairpin 
correlator that we measure is 

Hr(t) = ^(^T^TrLr^i + iOTrM^i)) (12) 

ti x y 

where the single fermion loop is 

L r (x,t)= $(xi-x)$(x 2 -x)(ip(x u t)rij(x2,t)). (13) 

x ly x 2 

r is a product of Dirac matrices: 75 and 7075 are studied for the eta prime. 
The connected correlator uses the same (separable) weighting function 

Cr(t) = (^2 E Hxi-xMx 2 -x)i;(x 1 ,t)T^(x 2 ,t) £ 0)r^(y 2 , 0)). (14) 

X 11,12 2/1,2/2 

We take the weighting function to be a Gaussian, <i>(x) = exp(— |a;| 2 /rg), of width ro = 3a. 
In coordinate space we fit 

C(t) = (exp(-mt) + cxp(-m(T - t)) (15) 

2m 

and, since 



4 



the coordinate space hairpin is 



Z u 2 

H(t) = — -j^((l + mi) exp(-mt) + (1 + m(T - t)) exp(-m(T - t))). (17) 

We break up the computation of the single loop into two terms. The first term includes the contribution of the low 
eigenmodes of D, and is done exactly (for those modes): Chiral eigenmodes of 13(0)^13(0) are labeled as 4>j±(x, t), the 
plus or minus sign corresponding to the chirality. They have associated eigenvalues Xj. These modes are recoupled 
into eigenmodes of Z3 _1 (m). The only mode mixing occurs between degenerate eigenmodes of Z3(0)^Z3(0). The single 
fermion loop becomes 

N 

Li(a?,t) = 2 [ Tr a j+ & j+ (x,t)r% j+ (x,t) + Tr o^-tfj.fotjrtfj-.fot) > (18) 

3 

where &j±(x,t) = J2 Xl ®( x ~ x i)4 > j±{ x ii't) is the convolution of the jth eigenmode with the smearing function <E>. 
Here, defining /i = m/(2xo), tj = Aj/(2xo), the ctj±, the eigenvalues of Z3 _1 (m), are 



1 A»(l-e?)±*eiyl-e? 
ai± = 2x~ e 2 + ^(l-e]) ' (19) 

For zero modes only one term in the brackets in Eq. (|l^) exists. 

The second part of the loop is computed using a stochastic estimator. We cast a vector of Gaussian random 
numbers \r(x, f)) on every site of the lattice, projected the modes used in L\ from it (so \r±) — |r) — \j)(j\ r ))i 
convoluted |n) with the smearing function, and used this vector as the source for a propagator. A final convolution 
of the propagator with the source vector produced a noisy estimator for the difference L(x, t) — L\(x, t). We averaged 
over twelve random sources per lattice. Each source was broken into two opposite chirality pieces, so a total of 24 
inversions, restricted to a single chirality sector each, were done per configuration. 

In retrospect, our simple stochastic estimator did not produce a useful signal. We first generated a 40-lattice 
data set, on which we found the complete hairpin correlator, from both low modes and stochastic estimator of high 
modes. For the pseudoscalar hairpin, all of the signal is contained in the low eigenmodes, which we captured exactly. 
Even so, we did not have a good enough signal from 40 lattices to analyze the hairpin graph, so we generated an 
additional 40 lattices. On these lattices we did an ordinary spectroscopy calculation for the connected diagram and 
computed the hairpin from the low modes. For our version of the overlap, the cost of finding the eigenmodes in 
addition to constructing ordinary propagators is neglegible because the eigenmodes can be used to precondition the 
Conjugate Gradient calculation of propagators. At 24 inversions per lattice, a stochastic evaluation of the higher 
modes' contribution is about twice as costly as an ordinary propagator calculation. 

We did not experiment with more sophisticated noise reduction methods |24[ . 



III. RESULTS 



It is an interesting exercise to ask how the contributions of different sets of eigenmodes contribute to the pseudoscalar 
hairpin correlator. We find that the low modes completely saturate the hairpin graph at all mass values studied. 
Because the plots are rather cluttered, we show results for only one mass value, am q = 0.04, in Fig. p|: The full 
correlator, including the low modes treated exactly and the high modes computed using the noisy estimator, is shown 
by octagons. Completely overlapping these points is the correlator built of the lowest 20 eigenmodes of D^D. 

The zero mode piece is shown by diamonds. At large quark mass there is a rather strong interference between the 
zero mode contribution and the nonzero modes, which pulls the full propagator down. Since the low modes themselves 
saturate the correlator, this interference must be dominated by them. 

Finally, the contribution of the nonzero modes by themselves is positive at short distance but becomes negative above 
a separation of t > 7 lattice spacings. We have displayed this by changing the plotting symbol for this contribution 
from crosses, when the signal is positive, to bursts, where it becomes negative. 

It is easy to understand why nonzero modes contribute negatively at large separation. The jth pair of nonchiral 
modes contribute a term to the hairpin of 

= ^EEE <w<^- t + t')n f (z',0) (20) 

V x.x' 
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where (see Eq. (|l8|)) 



Qj(x, i) = $(0;! 



(21) 



is a local smearing ol the chiral density of the jth mode. The nonzero modes have zero overall chirality and so 
at large distance the correlation function of the chiral density must become negative. This can be seen even on a 
single configuration. This behavior has been observed in the eigenmode studies of Ref. |2^|. It may be more than a 
coincidence that these authors saw the chirality autocorrelator (for the same data set as we studied here) becoming 
negative around a distance of 6 or 7 lattice spacings. 

Could we have used fewer modes? Fig. || shows the contributions to the hairpin correlator of the lowest 5, 10, 15, 
and 20 modes of H(0) 2 , for the am q — 0.04 data set. The pictures for the other masses are virtually identical. We 
could have kept only five modes of H(0) 2 and almost completely saturated the hairpin correlator, in our particular 
simulation volume. 
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FIG. 2. Contribution to the pseudoscalar hairpin for the am q = 0.04 data set from various subsets of eigenmodes: In all 
figures the full propagator is shown by octagons. In panel (a), the zero mode contribution is given by diamonds. In (b) the 
contribution of all nonzero modes is shown by crosses (where the contribution is positive) and bursts (where it is negative). In 
(c) the contribution of the lowest 20 modes of D is shown by squares. 
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FIG. 3. Saturation of the pseudoscalar hairpin graph at am q = 0.04 by 5, 10, 15, and 20 eigenmodes of H(0) 2 (with symbols 
octagons, diamonds, squares, and crosses). Bursts show the full correlator. 

We have performed correlated fits to the connected and hairpin pseudoscalar graphs to Eqs. (|l5|) and (|l7| ) and 
extracted lattice predictions for /j^/Nf. An example of a fit, to the am q = 0.04 low-eigenmode hairpin, is shown in 
Fig. ||. The error bars on the points are the naive ones and do not show the considerable correlation between the 
data from different time slices. A set of fits, for a time slice range of t m i n = 3 out to t rnax = 12, is shown in Fig. 
^. We have separated fits to using hairpin correlators with truncated quark propagators and fits using the full quark 
propagator in the hairpin. All of the fits to the low-mode hairpin have good confidence levels for t m in > 5 — 6. Pion 
masses are also a byproduct of the fit. Their values are stable for t m i n > 3 — 4 and equal within uncertainties to the 
results of fits to pion propagators in isolation. 

The "full" signal is much noisier than the low mode truncation (even more than one would expect, knowing that 
the low mode data set is twice as big), and, since both are equal within uncertainties, we quote the latter. Our chosen 
best fits, typically over a range 6-12, are shown in Table || and plotted in Fig. |^. We have also extrapolated our results 
for fiy/Nf to zero quark mass, assuming that fj^/Nf is a linear function of the quark mass, using a single elimination 
jackknife. 

Taking the lattice spacing from the rho mass, i.e. a = 0.13 fm, and setting Nf = 3 we predict /io = 770(54) MeV 
x(a~ 1 /1520 MeV). Had we used the Sommer parameter to set the scale, i.e. a — 0.11 fm or a -1 = 1770 MeV, the 
mass would have come out about 16% higher. This difference is, of course, a reflection of the inherent scale uncertainty 
in quenched QCD. Evaluating the left-hand side of Eq. (Q) with observed particle masses gives [1q = 860 MeV. 
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FIG. 4. Connected correlator (crosses) and hairpin (octagons), for the lowest 20 eigenmode hairpin, at bare quark mass 
am q = 0.04, showing the result of a correlated fit to the two correlators over the range t — 6 to 18 (folded to 6 — 12). 
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FIG. 5. Extracted pseudoscalar hairpin coupling Ho/Nf from fits from distance t to 12 at am q — 0.04. Octagons show the 
results of fits to the hairpin correlator when the quark propagator is truncated to the lowest 20 eigenmodes (80 lattice data 
set) of D^D, and crosses show fits using the full propagators (40 lattice data set) in the hairpins. 
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FIG. 6. Pscudoscalar hairpin coupling fio/Nf in lattice units as a function of quark mass. The jackknife-extrapolated 
zero-quark-mass value is shown as the burst. 
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FIG. 7. Contribution to the am q — 0.04 7075 hairpin from low modes of D: 5 modes (small diamond with cross when 
positive, plus sign when negative); 10 modes (diamond and fancy cross), 15 modes (square and burst) and 20 modes (octagon 
and cross). 

If the "standard picture" of the hairpin graph makes sense, we should be able to measure /j,q using any interpolating 
field at the source and sink points. In particular, the axial vector (7075 — 7075) hairpin correlator is an interesting 
operator to study. No zero modes contribute to its hairpin correlator, and yet the pseudoscalar and axial vector 
hairpins should both give the same result for /ijj - the only difference in the two channels is that the coupling of the 
external current source to the pseudoscalar channel is different. 

We attempted to measure /Zg from this correlator, without success. We raise the point as a potentially interesting 
exercise if future high-statistics data become available. 
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At short distances this correlator is negative. At larger distances it swings positive, before disappearing into noise 
at large separation. (See Fig. fijfor the am q — 0.04 data set.) The node in the correlator makes it impossible to fit the 
data to the simple form of the hairpin we have used above. Notice from Fig. fj] that the hairpin correlator saturates 
more slowly with number of eigenmodes, than the pseudoscalar correlator does: this, plus the nodes in the signal, 
makes the analysis impossible to do. 

Now we return to the pseudoscalar correlator. To make contact with other chiral observables we need a lattice 
determination of We do this by measuring the pseudoscalar density P (as one vertex of a two-point function) and 
extracting / ff from m^f n — 2m 9 (0|P|7r). Because the overlap action is chiral, the renormalization factors for the quark 
mass and the pseudoscalar density cancel. The pseudoscalar density is taken to be the naive operator P(x) = ipj^if). 
Note that because we use the "subtracted" propagator D^ 1 , Eq, (||), our measurement is equivalent to the use of the 
"order a 2 improved current" -0(1 + D/(2xo))7s(l + D/(2xo))ip. The results of this exercise are shown in Fig. || and 
also given in Table |[ Using again the lattice spacing from the rho mass we find in the chiral limit, obtained by linear 
extrapolation in the quark mass, f n = 163(21) MeV. Had we used the Sommer parameter to set the scale we would 
have obtained f„ = 193(25) MeV. The latter number is obviously much larger than the experimental value of 132 
MeV. We prefer the lattice spacing from the rho mass for phenomenology with the overlap fermions. A comparison 
of the point-to-point correlator in the pseudoscalar channel with instanton liquid model predictions also requires a 
large lattice spacing of 0.13 fm. This correlator is proportional to (0|P|7r) 2 . 
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TABLE I. Table of best-fit lattice parameters a 2 fio/Nf, af^, a 4 \, and 8 as a function of quark mass and after jackknife 
extrapolations to zero mass. 
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FIG. 8. Pseudoscalar decay constant (in lattice units) from this overlap action. Crosses are our data and a jackknife linear 
extrapolation to m q — is shown by the square. 



The two parameters which are related to $ are the zero mode susceptibility \ — Mo/w/C^-^/) an< ^ the P a ~ 
rameter which characterizes the strength of the hairpin contribution in quenched chiral perturbation theory Jl2[ , 
6 = I (8ir 2 N f f%) . We compute each of these parameters for each value of quark mass using a combination of 
single-elimination jackknife with correlated fits to pairs of two-point functions, and extra pol ate the results linearly in 
the quark mass to m q — 0. The resulting values are shown in Table | and in Figs. ^ and [l0[ 
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FIG. 9. The quenched chiral perturbation theory parameter 8 — fig / (Sir 2 N / /^) evaluated for Nf = 3. Crosses are our data 
and the extrapolated value at m q = is shown by the octagon. 
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FIG. 10. The inferred zero mode susceptibility in lattice units from the hairpin graph \ = / i oiW(4A'/) evaluated for Nf — 3. 
Crosses are our data and the extrapolated value at m q = is shown by the octagon. The susceptibility measured directly from 
zero modes is shown with the square. 

We can also perform a direct measurement of the topological susceptibility by counting zero modes: the topological 
charge Q on each configuration is just denned as the difference in number of opposite chirality zero modes in the 
configuration. For our 80 lattices, we find (Q) = 0.31 ± 0.32. We can (formally) eliminate the (Q) 2 term from 
the average by combining our data with a parity-reversed copy of every lattice. (This is equivalent to the usual 
replacement of a propagator by its real part in conventional spectroscopy or matrix element calculations.) The zero 
modes give us the result (Q 2 ) = 8.24 ± 1.11, and 
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a X 



V 



1.98(27) x 10" 



(22) 



This result agrees with the calculation of \ from /Xq. The quality of our signal is similar to what is seen in calculations 
of x using pure gauge observables (lattice analogs of the topological charge density) with similar statistics. Converting 
this lattice number to a continuum result using the Sommer scale, as is conventional for the topological susceptibility, 
and taking the fourth root yields x 1 / 4 = 213(7) MeV from the direct counting of zero modes and x 1/4 = 226(12) MeV 
from the hairpin fit. Pure gauge studies produce similar numbers to what we see [p6| . 

Obviously, with the scale set by the rho mass, the result would be about 16 per cent lower, 180(6) and 191(11) 
MeV. 

Our measurement of 8 ~ 0.093 is a bit higher than recent determinations made using nonchiral actions [p|, 28| . f27| . 
which measure S ~ 0.065 (but somewhat smaller than what one would infer from the W-V relation with physical 
masses, Eq. (Q)). The first calculation is done using clover fermions at (3 = 5.7. Clover fermions have real eigenmodes 
which make a nonvanishing contribution to Trj^D^ 1 . The overlap action counts each zero mode with unit weight. 
It is our experience that as the value of the real eigenvalue of a clover mode moves away from zero, the associated 
chirality also decreases, and so it is not surprising that our result is larger. The other two calculations determine S 
from the small quark mass behavior of the pion mass, i.e. from observing quenched chiral logarithms, and can easily 
have quite different systematic uncertainties than our measurement. 

We do not see any effects of quenched chiral logarithms in any observables. A fit to the pion mass in various 
channels (pseudoscalar, axial vector, difference of pseudoscalar and scalar) to the form suggested in Ref. p7] , m 2 . = 
Cm q (l — 5 log(Cm 9 /A 2 )), chooses 8's in the range 0.0-0.2, but with uncertainties from any fit also of the same order. 
The main effect of quenching we have seen is the presence of zero modes in the pseudoscalar and scalar channels; in 
contrast to the hairpin case, these are finite volume artifacts. 

We do not have any useful signals from other channels. The scalar channel is dominated by zero and low eigenmodes, 
but the major structure in the channel is just the constant chiral condensate. In the vector channel the low-eigenmodc 



12 



correlator even has the wrong sign-it is negative over most of its range (compare Fig. [0]). The absence of a signal 
could be consistent with a naive expectation from Zweig's rule, namely that the size of the hairpin in the vector 
channel is small. The low eigenmodes which we include exactly do not seem to make much of a contribution to 
disconnected diagrams in the vector channel, either. This feature is expected in instanton liquid models 
We would feel much more comfortable making these statements if we had a real signal p3] . 
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FIG. 11. Contribution to the am q = 0.04 -ji hairpin from the lowest 20 modes (square when positive, burst when negative) 
as compared to the full correlator (octagon when positive, cross when negative). 



IV. CONCLUSIONS 



For overlap actions, and in the quenched approximation, the Witten-Veneziano formula is an exact relation between 
the eta-prime mass inferred from the hairpin graph and the fermionic zero-mode susceptibility. However, the physics 
of the Witten-Veneziano relation in the overlap is rather different than in "standard derivations:" the eta-prime is not 
a real particle; the mass we measure is the size of a quenched-artifact coupling between two flavor-singlet Goldstone 
bosons. For quenched QCD, the gauge configurations which "fill in the white space" in the hairpin diagram are the 
ones which produce zero modes. We have not addressed the question of what choice of contact term is needed to 
equate a particular definition of the topological susceptibility to the zero mode susceptibility. 

The numbers we have found for the inferred eta-prime mass and the zero mode susceptibility are not all that 
different from previous results using nonchiral fcrmion actions. We believe that the theoretical underpinning of our 
calculation done with a chiral action is more reliable than any calculation done with a nonchiral action. We are well 
aware that our calculation is performed at only one value of the lattice spacing, and that simulations at several lattice 
spacings are necessary for an honest extrapolation to a continuum value. 

The strong coupling of low eigenmodes to the hairpin amplitude allowed us to perform the numerical simulation. 
This connection between the zero mode susceptibility and the coupling strength is not automatic for a nonchiral 
fermion action. Nevertheless, we expect that small eigenmodes of the Dirac operator will make a large contribution 
to the pseudoscalar hairpin correlator. We suggest that future studies of the hairpin graph, even done using nonchiral 
actions, can reduce the measurement noise by first finding all the zero or near-zero eigenmodes of the Dirac operator 
so that their contribution can be included exactly |3l[| . 
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